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ABSTRACT 

The recent discovery of persistent GRB counterparts at lower frequencies permits several 
important conclusions to be drawn. The spectrum of GRB970508 is not consistent with an 
external shock origin for both the prompt GRB and the persistent emission, suggesting that 
at least the prompt radiation is produced by internal shocks. Comparisons among three GRB 
with counterparts (or upper limits on them) establishes that GRB are not all scaled versions 
of similar events. The angular size inferred from the apparent observation of self-absorption 
in the radio spectrum of GRB970508 a week later implies that its expansion had slowed to 
semi-relativistic speeds. This permits a remarkably low upper bound to be placed on its residual 
energy, suggesting either that radiation has been > 99.7% efficient or that the initial outflow was 
strongly collimated. Observations of self-absorbed radio emission from future GRB may permit 
direct measurement of their expansion and determination of their parameters and energetics. 
We estimate initial Lorentz factors for GRB970228 and GRB970508 70 ~ 100, and present 
a solution for the evolution of a blast wave with instantaneous cooling. The visible flux is 
predicted to rise initially as a power of t slightly less than unity, in agreement with observations 
of GRB970508. 


Subject headings: Gamma Rays: Bursts — Gamma Rays: Theory 


1. Introduction 

The recent discovery by BeppoSAX and coordinated visible, infrared and radio observations of 
persistent emission from gamma-ray bursts (GRB) has answered several important questions. Accurate 
coordinates of visible counterparts have led to conclusive evidence that GRB are at cosmological distances, 
definitively resolving the oldest and most fundamental question in GRB astronomy. Yet other questions, 
such as the nature of the GRB emission mechanism, remain unanswered, and new questions have arisen. 
Most important, it is unclear how the persistent emission is produced and what is its relation to the prompt 
soft gamma-ray radiation which has defined GRB since their discovery in 1973. 

In this paper we examine two different classes of models for the persistent emission. In one class it 
reflects continuing gamma-ray activity, not directly observed because of the comparatively low sensitivity 
of GRB detectors to emission of low intensity but long duration. In the other class the persistent emission 
is the long-predicted afterglow produced when a relativistically moving shell sweeps up ambient matter and 
its large initial Lorentz factor decays as its energy is shared and radiated. We also examine two different 
classes of models for the prompt GRB emission. In one it is produced by “internal shocks”—interaction 


1 Permanent address: Department of Physics and McDonnell Center for the Space Sciences, Washington University, St. 
Louis, Mo. 63130 



- 2 - 


between different shells each moving at different but relativistic speeds. In the other class it is the result 
of an “external shock” produced when a relativistic shell collides with the ambient matter. We will argue 
that the presently available observations (chiefly of GRB970111, GRB970228 and GRB970508) present 
strong evidence that the prompt GRB emission is that of internal shocks. The explanation of the persistent 
emission is less clear, with some evidence for continuing gamma-ray activity and some for afterglows; it 
is possible that both processes contribute with differing importance in different GRB. We also present 
evidence that in at least one burst (GRB970508) either the external shock is an extremely efficient radiator 
of the kinetic energy or the initial mass outflow and gamma-ray emission are strongly collimated. 


In §2 we summarize the two classes of models for persistent emission. If the prompt GRB were 
produced by an external shock then in a burst whose persistent emission is afterglow it should be possible 
to extrapolate the properties of the afterglow back to the prompt GRB. In §3 we present evidence that this 
cannot be done for GRB970508, implying that its prompt emission was the result of internal shocks. This 
conclusion is particularly interesting because the prompt emission of GRB970508 had a simple single-peaked 
time profile, consistent with an external shock; external shocks are not satisfactory sources of the gamma-ray 
emission even where permitted by the time-history. We compare three bursts in §4 and conclude that 
GRB are not all scaled versions of similar events. In §5 we discuss evidence for self-absorption at radio 
frequencies in GRB970508, and estimate the size of the emitting region. From this we can infer either that 
radiation has been very efficient or that the initial expulsion of relativistic matter was strongly collimated. 
Collimation differing from burst to burst is one possible explanation of the comparisons presented in §4. 

In §6 we interpret the pulse shape and spectral history of GRB970228 and GRB970508 as the results of 
internal shocks followed by external shocks, and estimate their parameters. §7 presents simple estimates 
and scaling laws for the properties of afterglows produced when radiation is efficient, as suggested here. §8 
contains a summary discussion and considers the prospects for determining the parameters of GRB and 
their afterglows by using measurements of self-absorbed synchrotron radiation to measure their expansion. 
In an appendix we rederive and explain the result for the self-absorbed intensity of a GRB given (with little 
explanation) by Katz 1994a. 


2. What Makes Persistent Emission? 

Gamma-ray bursts are generally agreed to be associated with the expulsion of relativistic debris by a 
condensed object. The fundamental argument for this is the presence of ganrma-rays with energies greater 
than m e c 2 . If gamma-ray bursts are distant and luminous (the evidence for cosmological distances and 
great luminosities is now compelling) such ganrma-rays can avoid destruction by 7-7 pair production only 
if they are narrowly collimated in a radially outward direction and their radius of emission is large enough. 
A sufficient, and probably necessary, condition for this to occur is that they be emitted by matter moving 
radially outward at relativistic speed with a Lorentz factor 1. 


2.1. Afterglow Models 

The outward moving matter will slow as a result of its interaction (generally assumed to be collisionless, 
for the collisional mean free path is extremely long) with any surrounding medium. Shocked matter will 
radiate, and as the Lorentz factor of the debris decreases the radiation will shift to progressively lower 
frequencies and the corresponding time scales will lengthen. 
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Considerations of this kind led several authors ( [Katz 1994a , Katz 1994b, Meszaros & Rees 1997, [gari 
fe Piran 1997a , Vietri 1997 , Waxman 1997 ) to predict that GRB would be followed by an afterglow at 
frequencies gradually declining from soft gamma-rays through X-rays to visible light and radio waves. 

It is possible to make many different assumptions concerning the efficiency of radiation, reconversion of 
post-shock thermal energy to bulk kinetic energy, collimation of the initial debris and development of the 
magnetic field. Meszaros & Rees 1993 and Meszaros, Rees & Papathanassiou 1994, who integrated GRB 
spectra through the stages which are now called afterglow, and Paczynski fe Rhoads 1993| , who considered 
radio emission at very late times, demonstrate the possible variety and complexity of these assumptions. 
As a result, many different predictions of the quantitative development of the radiation spectrum with time 
are possible. 


All such models in their basic form, in which only one set of assumptions is made, have the general 
property that characteristic time scales and peak intensities are power law functions of frequency. The 
reason for this is that if there is no characteristic Lorentz factor defined between the Lorentz factor 70 of 
motion initially produced by the source and unity, then there is no opportunity for a break in the power 
law dependence of any function on any variable between the values defined at these limits; such a break 
would define a new characteristic Lorentz factor, which by assumption is not present. Of course, it is 
possible to define more complex models in which this assumption does not hold (for example, radiation may 
be assumed efficient in some range of Lorentz factors, and inefficient outside that range). This possibility 
opens the door to a plethora of parameters and models of unmanageable complexity. Such models may fit 
the data, but it is hard to decide if this is because they are physically correct or because it is easy to fit 
multiparameter functions of no physical significance; they will not be considered further. 


Afterglow models generally predict that the intensity at a given frequency should rise until a 
characteristic time, and then decay (Katz 1994a, Katz 1994b|). This appears to have been observed 


for visible light from GRB970508 (Castro-Tirado et al. 1997, Djorgovski et al. 1997 , Galarna et al 


1997b, Jaunsen et al. 1997, Schaefer et al. 1997), although the initial rise is defined by only a very few 


observations. The characteristic time is a decreasing function of frequency, short for gamma-rays, longer 
for visible light and yet longer at infrared and radio wavelengths. Before this single maximum the intensity 
should rise monotonically, and afterwards it should decay monotonically, without fluctuation, as a smooth 
(usually power law, when well away from the maximum) function of time. This is so far consistent with 
nearly all the data (see, however, Kopylov et al. 1997 ) but they are so sparse that this conclusion can 
hardly be said to have been proven. 


2.2. Continuing Emission Models 


In an alternative class of models ( Katz, Piran fe Sari 1997 ) the persistent radiation at lower frequencies 
(X-rays, visible light, infrared and radio waves) is produced by the same processes which make the initial 
GRB. The gamma-ray emission continues at lower intensity, perhaps maintaining a constant ratio to the 
lower frequency intensity. 


There is little direct evidence for or against models of this class. Emission at frequencies below X-rays 
has never been observed simultaneously with a GRB for instrumental reasons; only high upper bounds can 
be set on its intensity. Analogously, weak continuing gamma-ray emission is nearly unobservable by a GRB 
detector, which is designed to trigger on intense short bursts. Continuing emission models agree with the 
spectral slope of persistent emission in I, V and X-ray bands in GRB970228 (|Katz, Piran fe Sari 1997). 
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Some indirect evidence suggests that gamma-ray activity may continue long past the nominal duration 


(usually < 100 s) of GRB. High energy gamma-rays were observed 5000 s after GRB940217 (Hurley et al 


1994). The “Gang of Four” spatially coincident (to within experimental accuracy) bursts of October 27-29, 


1996 (|Meegan et al. 1996 , Connaughton et al. 1997|) indicate repetitive activity from a single source, which 
can equivalently be described as a single burst lasting two days, with brief periods of intense emission 
amidst a longer period of much weaker or undetectably faint emission. 


Models of this class predict that continuing gamma-ray activity should accompany the lower frequency 
radiation. In analogy to the behavior of the intensity during observed shorter GRB, the intensity should 
fluctuate irregularly in all bands, although this may be hard to detect in faint objects for which observations 
require extended integrations. 


3. Evidence for Internal Shocks in GRB970508 


We wish to test the applicability of the models discussed in §2 to GRB. For a variety of reasons (satellite 
pointing constraints, instrumental sensitivity, Earth occultation, daylight, moonlight, etc.) observations of 
persistent emission by GRB are sparse. In order to test continuing emission models it is necessary to have 
data obtained simultaneously at several frequencies. The spectrum of GRB970508 does not support these 
models. 

GRB970508 is therefore an excellent candidate for afterglow models. All the soft X-ray data and nearly 
all the visible data were obtained after the peak intensities in those bands, and therefore reflect the behavior 


of a high energy tail to the particle distribution function ( Wijers, Rees fc Meszaros 1997 ) and not the 
evolution of its characteristic energy or gross energetics. Fortunately, if there are sufficient data to define 
the peak flux F maxu at a frequency v this quantity can be determined without the need for simultaneous 
measurements across the spectrum, or for frequent sampling on both sides of the maximum. Therefore, it 
may be easier to compare F maxv with theories than the instantaneous F v or than F v (t) at a single u. 

The data available for GRB970508 are shown in Figure 1. The radio data are shown as a lower 


limit because (in the most recent available data; Frail et al. 1997b) the flux is still rising. 

The visible 

(Castro-Tirado et al. 1997, Djorgovski et al. 1997 

, Galama et al. 1997b, launsen et al. 1997 

, Schaefer et al. 

1997), hard X-ray (Costa et al. 1997b) and gamma-ray (Kouveliotou et al. 1997) maxima were observed. 

The soft X-ray flux (Piro et al. 1997) is a 

lower limit because its maximum was missed. Some specific 

afterglow models (Katz 1994a, Katz 1994b, 

Meszaros & Rees 1997) predict F maxv to be independent of v; 


these are evidently excluded. 

These points are not consistent with a single power law. Simple afterglow models (§2.1) cannot describe 
the entire spectrum, regardless of the specific assumptions made. It is possible to resolve this problem 
in many ways, for example by introducing two different regimes into the afterglow model, with suitably 
chosen parameters. We suggest, instead, that the brief measured hard X-ray and gamma-ray emission are 
not produced by the same process as the persistent emission at lower frequencies. |5ari fe Piran 1997a| 
argued that the gamma-ray emission of multi-peaked GRB could not be the result of an external shock, but 
rather was the result of internal shocks. We suggest that this conclusion is also applicable to GRB970508, 
resolving the discrepancy with afterglow models, which may explain the lower frequency emission. Because 


GRB970508 was single-peaked in gamma-ray intensity (Kouveliotou et al. 1997) the arguments of Sari fc 


Piran 1997a do not apply to it directly. Our argument for internal shocks in GRB970508 then suggests that 
in all GRB gamma-ray emission is produced by internal shocks, whatever the shape of their time profiles. 
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4. GRB Differ Qualitatively 


It is natural to ask if all GRB are scaled versions of similar events. External shock models predict 
that they are, because external shocks may be scaled by powers of their basic parameters: distance, total 
energy, initial Lorentz factor and ambient gas density. Table 1 presents ratios of F maxl/ in four bands 
between GRB970508 and GRB970228. These two bursts have quite different spectral properties. A similar 
conclusion may be reached by comparing GRB970508 to GRB970111 (Table 2). The different spectral 
properties are again evident. 


These tables summarize the obvious fact that GRB970508 was a comparatively weak burst in 
gamma-rays. It was brighter in V than GRB970228 and brighter at 1.43 GHz than GRB970111, two bursts 
which were much stronger in gamma-rays. The obvious explanation is the same as that reached in §3: the 
gamma-ray emission is the result of a different process than that which produces the persistent emission at 
lower frequencies. 


The weakness of GRB970508 in gamma rays poses other problems. Its gamma-ray fluence (|Kouveliotou 
et al. 1997 ) and redshift ( Metzger et al. 1997 ) imply an energy output (if isotropic) of ~ 3 x 10 51 erg. 


This makes it difficult to argue that its comparatively low gamma-ray fluence results from unusually 
inefficient gamma-ray emission, because if it had its observed V magnitude and the same V/gamma-ray 
ratio as GRB970228 the observed fluence would imply (assuming isotropy) the radiation of ~ 10 53 erg of 
gamma-rays. This is excessive; it is unlikely ( Katz 1997 , Kobayashi, Piran fe Sari 1997| ) that internal shocks 
are even 50% efficient! 


One possible explanation is that GRB970508 is, instead, an unusually efficient source of visible and 
radio emission. This would require that the efficiency of lower-frequency afterglows varies greatly from 
burst to burst, which seems implausible; a priori , one might expect more variation in the efficiency of a 
complex process like internal shocks, which depend on the detailed temporal dependence of the Lorentz 
factor, than in that of a simple one like a blast wave, which only depends on a few parameters. 

A more promising alternative is that the gamma-ray emission of GRB970508 was beamed towards us, 
so that the total power radiated was much less than that implied by the assumption of isotropic radiation. 
Beaming could solve all energetic problems, but the flux ratios require that other bursts show an even 
greater degree of gamma-ray beaming in order to explain their comparatively (to GRB970508) greater 
strength in gamma-rays in comparison to visible or radio fluxes. Further evidence for beaming is discussed 
in the next section. 


A third possible explanation is that discussed in §2.2. If the lower frequency persistent emission 
(from soft X-rays down to radio waves) is the product of continuing activity by the compact source of the 
gamma-ray burst (Katz, Piran & Sari 1997), then it would not be surprising if the strength of this activity 
differed from burst to burst, just as the time dependence of the gamma-ray emission differs. 


5. Self-Absorption, Radiative Efficiency and Beaming in GRB970508 


The observed (Frail et al. 1997b) radio spectrum of the persistent emission of GRB970508 on May 
15, 1997 was oc v 11 . This is an extraordinarily steep slope, and is naturally interpreted as the transition 
between optically thick and thin conditions in the frequency range observed. The expected self-absorbed 
spectrum ( Katz 1994a ) is oc v 2 (differing from the familiar i/ 5 / 2 of self-absorbed synchrotron radiation 
because in relativistic shocks the nearly all the low frequency radiation is emitted by electrons with a single, 
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and much higher, characteristic synchrotron frequency). The expected spectrum at frequencies between the 
self-absorption frequency and the characteristic frequency is oc z/ 1 / 3 . 


A theory for self-absorption in GRB was presented by Katz 199~4a| (pp. 255-256) and is given in more 
detail in the Appendix. The result is that below the self-absorption frequency 


F v 


it 

2Trv 2 m p ((l + z)jp, 


(1) 


where R is the radius of the radiating shell, 0 its redshift and D = (2c/H 0 )[l — (1 + z) -1 / 2 ] its proper 
distance. Here we include the factor £ = ksTe/("fm v c 2 ) describing the degree of electron equipartition in 
the plasma shock-heated to an internal energy per particle ym^c 2 and moving with Lorentz factor 7 . Using 


the measured (Frail et al. 1997b) flux density on May 15, 1997 at 1.43 GHz as the optically thick value 
(this data point is only significant at the 2 er level, but is consistent with the very significant extrapolation 


of higher frequency flux densities) and taking H 0 = 60 km/sMpc and 2 = 0.835 ( Metzger et al. 1997 ) 
as the actual redshift yields R sa 4.1 x 10 16 (/ _1 / 2 cm « 1.2 x 10 17 cm for £ = 1/9 (a likely upper limit, 
assuming equal energy in electrons, ions and magnetic field). The usual result for the size after an interval 
t is R ~ 2y 2 cf, which would yield 7 Ra 2, for which the relativistic approximation is barely valid. Sari 1997 
pointed out that because of the gradual deceleration of the shell the correct expression is R r\j 8 y 2 ct (if 


7 oc r , as in the blast wave model of Katz 1994a and in §7, this becomes R ~ 147 2 ct), implying that the 


expansion is nonrelativistic by the time of these observations, a week after the burst. 

This may be hard to reconcile with the inferred minimum energy ~ 10 52 erg because at a density 
of n = 1 cm -3 only 1.5 x 10 28 g (1.4 x 10 49 erg of rest mass energy) is contained within a sphere of 
R = 1.3 x 10 1 ' cm. In order for the motion to have been only marginally relativistic the kinetic energy 
must have been reduced to ~ 10 49 ni erg from its initial value, where ri\ = n/(lcm -3 ), requiring a radiative 
efficiency e > 99.7%, and probably at least 99.9%. Such efficient radiation is not possible from internal 


shocks for kinematic reasons (Kobayashi, Piran fe Sari 1997, Katz 1997), but may be in an external shock 


( Katz 1994a , Sari, Narayan fe Piran 1996 , Sari 1997 ). The hypothesis of extremely efficient radiation may 
be tested if a complete energy budget is obtained for future GRB, including the early stages of the external 
shock and the difficult spectral range of soft X-rays and the far ultraviolet; the early stages of the afterglow 
must emit as much radiation as the GRB itself, and perhaps much more. 


Alternatively, if the relativistic debris shell is initially only produced within a narrow solid angle Q 
its energy requirement is reduced. Once its motion becomes nonrelativistic the collimation is largely lost, 
as shock-heated matter expands in all directions. Then the inferred ~ 10 49 ni erg of kinetic energy may 
represent most of the actual burst energy, not requiring efficient radiative loss. The measured GRB fluence, 
now taken into the solid angle H, implies a radiated energy of ~ 2.5 x 10 5 °H erg. Then 

H « 0.04—-—m. (2) 

1 — e 


In internal shock models plausible values of e 


0.5 imply a collimation opening angle 9 


0.1 n[ 


1/2 


radian. 


_1 / Q 

This degree of collimation requires 70 )> 1/9 s=s lOrij , consistent with other estimates. It also argues 
against n\ <C 10 -3 , such as would be found in galactic halos or intergalactic space, for this would lead 
to excessively small 9 ; even n\ = 10^ 3 , as found in much of the volume of the Galactic disc would imply 
7 > 300. 
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6. Possible Parameters of GRB970228 and GRB970508 


6.1. GRB970228 


The gamma-ray (40-600 keV) light curve of GRB970228 displays a strong and rapidly variable peak 
about 4 seconds long, followed by a much weaker and smoother peak that begins about 40 seconds later 


( Costa et al. 1997a ). The first pulse is variable on a time scale of ~ 0.1 s ( Hurley et al. 1997 ). The 
X-ray (2-30keV) light curve of GRB970228 displays two comparable peaks, corresponding in time to the 
gamma-ray peaks. The observed flux in the second X-ray peak is comparable to the flux obtained from 
backward extrapolation of the persistent X-ray emission light curve. 

This suggests a simple interpretation: The first GRB pulse was produced by internal shocks while the 
second pulse signals the beginning of an afterglow produced by an external shock when the debris sweeps 


up ambient gas. Behavior of this kind was suggested by Sari fc Piran 1997a who pointed out that internal 
shocks could not convert the total kinetic energy to radiation and that a significant fraction of the kinetic 
energy should be emitted later in other parts of the spectrum. This interpretation allows us to estimate 
some of the parameters of the event. 

The overall duration of the first pulse corresponds to the duration of intense activity of the “inner 
engine” that emitted the ejecta. The variability on a time scale of ~ 0.1 s corresponds to variability in the 
“inner engine” and sets an (unsurprising) direct upper bound on its size of ~ 3 x 10 9 cm. 

The second pulse corresponds to the the onset of significant radiation by the external shock. This 
takes place when the ejecta (with energy E, Lorentz factor and proper mass M = E/j qc 2 ) sweeps 
up an ambient mass M /70 at r = i ? 7 = [?>E / (AimmpC 2 ^)} 1 / 3 ~ 1.2 x lO 18 ^ 2 ^ 3 E^n^ 1 ^ cm, where 
1?52 = -E/(10 52 erg). The corresponding time of arrival of photons from the onset of the external 
shock is t e ~ Ry/(2 c7q) ~ 2 x 10 7 7 q 8 ^El ) !^n 1 1 ^ 3 s. Comparison with the time delay of 40 s yields 
140(l?52/ni) 1 / 8 . This is a reasonable value, but the dependence on E and n is so weak that no 


7o 

significant constraints can be placed on them. |Fenimore, Epstein fc Ho 1993[ |Woods fc Loeb 1995|, [Piran 


1996 give lower limits to 7 > 100 but Bari fc Piran 1997b , Katz 1997 give also an upper limit to 7 of the 
order of 10 3 for the observation of internal shocks. 


6.2. GRB970508 


The observation of a visible maximum 1-2 days after GRB970508 ( 

Castro-Tirado et al. 1997, 

Djorgovski 

et al. 1997, 

Galama et al. 19971, 

Jaunsen et al. 1997. 

Schaefer et al. 1997 

) suggests that at this time the 


characteristic synchrotron frequency of the afterglow passed through the visible spectrum. It should be 
possible to extrapolate this frequency backwards in time to t = t e , when 7 s=s 70 and the afterglow begins. 


In principle, this earlier part of the afterglow could have been observed by continuously monitoring 
GRB detectors. The fact that it was not seen in GRB970508 (in contrast to the suggestion made above for 
GRB970228) suggests that its peak frequency was too low for observation by GRB detectors, although we 
cannot be sure (without quantitative estimates of its radiative efficiency) how intense it would have been. 


Interpreting the nondetection by GRB detectors to mean that the peak characteristic afterglow photon 
frequency v < ^ 0 , we find 70 <) 60</ e ^ Cb 1 ^ 8 n 1 1 ^ 8 (/ir'o/100KeV) 1 / 4 . Although this may be less than lower 
limits estimated by consideration of gamma-gamma pair production ( |Fenimore, Epstein fc Ho 1993 , Woods 
fc Loeb 1995|, Piran 1996), these limits are smaller for afterglow at larger radii than for prompt GRB 











































emission, and are not directly applicable to a burst (such as GRB970508) for which no data exist on the 
high energy spectrum. It is also plausible that £ e , Cb or n\ are significantly less than unity, raising the 
upper bound on 70 . 


At the time of the visible maximum we can estimate 7^4 (Equation 7). Extrapolation (using 
7 oc f -3 / 7 ; §7) implies the blast wave will become only semirelativistic (7 < 2 ) after about a week, in 
agreement with the inference made in §5. 


The implied (using scaling laws from §7) X-ray (hv ~ 10 KeV) maximum occurred ~ 10 3 s after the 
burst, before the BeppoSAX follow-up observations (Piro et al. 1997) but well after the gamma-ray burst 
itself (Kouveliotou et al. 1997). For this reason the X-ray emission measured simultaneously (Costa ei 


al. 1997b) with the GRB was probably the low frequency extrapolation of the gamma-rays, produced by 


internal shocks, and not the beginning of the afterglow. The simultaneous X-ray emission also had F max 
much greater than the visible F max , inconsistent with the theory for both weakly (Katz 1994a, Katz 1994b| ) 
and strongly (§7) radiating blast waves, but consistent with the downward (oc iW 3 ; Katz 1994a, Kat? 


1994b, Cohen et al. 1997) extrapolation of the gamma-ray spectrum. 


7. Blast Waves with Instantaneous Cooling 


Previous relativistic blast wave models (Meszaros & Rees 1993, Katz 1994a, 

Katz 1994b, 

Meszaros, 

Rees & Papathanassiou 1994, 

Meszaros & Rees 1997, Vietri 1997, Waxman 1997) have generally assumed 


that only a small fraction of the energy is radiated, as in the Taylor-Sedov nonrelativistic blast wave theory. 
The arguments of §5 suggest that this may not be the case, at least in GRB970508, but that radiation may 
be very efficient, as indicated by the estimated short electron radiation times ( Katz 1994a , Sari, Narayan fe 
Piran 1996 ). Alternatively, the angular spreading of an initial narrowly collimated outflow as its Lorentz 
factor degrades will reduce the energy per solid angle with increasing distance along the axis of the flow in 
a manner qualitatively resembling the effects of efficient radiation. 


Here we develop a simple model for the evolution of a relativistic blast wave and its radiated spectrum, 
assuming that all the energy produced by the shock is instantly radiated. This is the limiting case opposite 
to that of negligible radiation, previously assumed. Strong coupling between the protonic and the electronic 
thermal energies is required. 


We describe the interaction of a relativistic blast wave with the ambient matter by a series of inelastic 
collisions between the debris and previously swept-up mass with proper mass M + m and Lorentz factor y(r) 
and infinitesimal shells of proper mass dm. The internal energy produced in each collision is instantaneously 
radiated. Conservation of energy and momentum yielcly : 


dj dm{r) 

7 2 — 1 M + m(r) ’ 

which can be integrated: 

7 + 1 = 7o + l / m(r) \ 2 

7 - 1 7o - 1 \ M ) 


(3) 

(4) 


2 Equations (71) and (72) of Katz (1994a) are wrong. 
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Then 


( 7 o + l)[(r/T) 3 + l] 2 + (7o-l) 
( 7 o + l)[(r/ J L) 3 + l] 2 -(7o-l)’ 


(5) 


where L = [ 3 M/( 47 rm p n )] 1 / 3 = [3E/(Anm p c 2 n 7 o)] 1/3 . L corresponds to the radius within which the ambient 
mass equals the proper mass of the debris. L is smaller than the Taylor-Sedov radius Rts within which the 
ambient rest-mass energy equals E and it is larger than i? 7 ; the ratios Rts '■ L '■ R~i = 7 o ^ 3 : 7(/ 3 • 1- 


There are three limits for this solution: 


7(r) 


7 o for r <C i? 7 

(■ r/L )~ 3 for 1?, < r < L 

l + 2(r/L ) 6 for L <C r 


( 6 ) 


Comparison of this solution with the Lorentz factor when there is no significant cooling reveals, naturally, 
that this Lorentz factor is smaller. For our purposes the intermediate range is the most interesting and we 
will use it in the following discussion. The 7 oc r -3 dependence found here is also that found if it is assumed 


the blast energy is distributed uniformly, or in a way which scales with r, as r increases (Katz 1994a), and 
a similar assumption is made in the Taylor-Sedov nonrelativistic blast wave theory. 

This Lorentz factor is also the typical “thermal” Lorentz factor of the relativistic shocked protons. 
The thermal Lorentz factor of the relativistic electrons is, by definition of the equipartition factor £ e , 

7 e = Ce( TO p/ TO e)7- The characteristic synchrotron radiation frequency from the forward shock (the 


backward shock crosses the shell and is gone before the shell enters the 7 <C 70 regime) is ( Piran 1994 ) 


hv : 


A A1 ~\T 4 1/2 

0.01 eV CeCe 7 W 


0.01eVC 2 C b / 2 ( y ) 


1/2 


(7) 


where the magnetic equipartition factor (b is the ratio between the magnetic energy density and the 
“thermal” energy density of the shocked protons. For typical parameters the characteristic frequency is in 
the visible range when 7 « 4, and the blast wave ceases to be relativistic when this frequency is in the 
infrared. 


The energy emitted at radius r and Lorentz factor 7(7-) will be observed around a time: 


t(r) = 


r/c 


14y 2 (r) 


= (L/14c)(r / L) 7 . 


( 8 ) 


Here we have replaced the usual factor of 2 in the denominator by 14 to allow for deceleration ( Sari 1997 ; 
see §5). 


Combining these results and setting the equipartition factors to unity we find: 


hv « 0.00011 eV(ct/L)~ 1 ' 2/7 n 1 1 /2 « 0.14 eV 


1 \ 12/7 £.4/7 / 7o V 

1 day ) 52 V300J 


■4/7 


-1/14 


(9) 


This gives a scaling law t oc v 'Z 12 which may be used to compare or predict the times of peak flux at 
different frequencies. Note also that 7 oc f -3 / 7 and r oc f 1 / 7 . 

We can now derive the scaling of F maxv and F v (t), using Equations ( 6 )-( 8 ): The time integrated 
fluence per unit frequency dE/dv = (dE/dr)(dr/dv) oc r 9 . This is emitted over a characteristic time oc r 7 , 
yielding 

Fmaxv oc T 2 (X V~ 1/6 . (10) 
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This should be valid throughout the regime in which the shell is moving relativistically, roughly 
corresponding to v from the infrared to X-rays. This law can, in principle, be compared with data like those 
shown in Figure 1, but in that Figure there is only one significant point in the applicable frequency range. 


At lower frequencies the instantaneous spectrum F v oc v 


1 / 3 (Katz 1994b), so that at a fixed frequency 


F v {t) oc F„ 


1/3 


i/3^6/7. 


( 11 ) 


The exponent of t is very close to the value of 4/5 obtained in a model in which radiation is inefficient 
Katz 1994^ , K qfal9946 . Observations of the persistent visible counterparts to GRB970228 ( |Guarnieri ei 
al. 1997 ) and GRB970508 ( Djorgovski et al. 1997 ) confirm the predicted rise before the maximum, but 
are not sufficient to determine the exponent quantitatively. The data for GRB970508 suggest an exponent 
of 0.7 ± 0.1, approximately consistent with the predictions, but this value may underestimate the true 
exponent because it is based on only two points, the second of which was obtained close to the maximum 
when the initial rise may have been leveling off. 


8. Discussion 


The persistent counterparts of GRB may offer as much complexity and variety as the classical GRB 
themselves. Future observations will provide a catalogue of their behavior, and will show if this is true, or 
if they may be described by a simple general model. Such a general description, with a few parameters 
to describe the possible collimation and degree of energy loss, as well and E, n and 70 , would support 
afterglow models, while complex “no two of them alike” behavior would support the suggestion of continuing 
emission. Afterglows are almost inevitable, and were predicted by many authors, but there is also evidence 
for continuing emission. Perhaps each may be important in different bursts. 


The apparent discovery of self-absorption in GRB970508 and the measurement of its redshift offer the 
prospect of detailed reconstruction of the history of a GRB blast wave. It may be possible to measure 
directly the blast wave radius as a function of time, and therefore to determine its slowing history. If 
radiation is efficient then by constructing an energy budget from the observed radiation it would be possible 
to determine E directly. By comparing the radiated power to the slowing history of the blast wave it would 
then be possible to determine n. It may also be possible to determine 70 (§6,7). 


If GRB are collimated, as suggested for GRB970508 in §5 and has been widely speculated on the 
grounds that their likely sources, coalescing compact binaries, are far from spherically symmetric, the 
analysis of afterglows will be much more complicated. They will have several more parameters, and will 
probably require numerical simulation. However, even the first clear demonstration of collimation would be 
significant. 


In Equation ( 8 ) most of the sensitivity of t to r comes from the 7 2 denominator; to a fair approximation 
t oc 7 ^ 2 . Combining this with v oc 7 4 (Equation 7) implies that t, will generally be proportional to 
approximately the -1/2 power of v. We expect this to hold for all fireball afterglow models, whatever 


their detailed physics. For example, this paper found t oc v 7 / 12 and the very different model of Katz 


1994a , Katz 1994b| found t oc z/ -5 / 12 . If 7 oc r -3 / 2 , as assumed in many models, then t oc zz -2 / 3 . In simple 
synchrotron radiative cooling models without hydrodynamics t oc . For similar reasons the exponent in 
Equation (11) is generally slightly less than unity. This suggests that it may not be possible to distinguish 
among fireball afterglow models on the basis of the time dependence of the characteristic frequency or the 
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pre-maximum rate of rise. However, these are not the only means available; for example, an energy budget 
will decide if radiation is energetically important. 


We thank Lifan Wang and Re’em Sari for discussions, Washington University for the grant of sabbatical 
leave, the Hebrew University for hospitality and a Forchheimer Fellowship, and the U. S.-Israel BSF, NASA 
NAG 52682 and NAG5-1904 and NSF AST 94-16904 for support. 


A. Self-Absorption 

A quantitative calculation of the self-absorbed flux from a rapidly moving emitter is not trivial, because 
a retarded-time calculation must be performed to allow for the fact that different points on the radiating 
surface are observed at different radii and retarded times. Deceleration and limb darkening introduce 
further complications. Here we only demonstrate that for shock-heated emitters the Lorentz factor does not 
appear in the expression for the observed flux density. 

In the co-moving frame of the shocked radiating matter moving with Lorentz factor 7 1 the 

temperature is given by T = (fm p c 2 /k B , where the Doppler shift factor / = (1 + vcosO/c )y ~ 2y ( Katz 


1994a took £ = 1 and used a variable 7 p which is here equal to / 2 ). In that frame the Rayleigh-Jeans 


brightness 


2jA 

r-2 


B„ = B—Qf mpC 2 . 


(Al) 


The observer measures a blue-shifted temperature C,f 2 m p c 2 /k B , and hence observes a brightness 

K = 2zz 2 C fm p . 


(A2) 


However, the observer sees significant radiation only from a solid angle « nr 2 /(f 2 D 2 ) because of the 
relativistic beaming of radiation into a cone of angle 1//. The observed intensity 


F v 


2Trv-Cm p j^- 1 


(A3) 


the powers of / cancel. 


An alternative derivation begins by noting that the quantum occupation number n v is a Lorentz- 
invariant scalar. The intensity 

I v = v 3 n v . (A4) 


When Doppler shifting to the observer’s frame /„ 
number 


■ Lf 3 - 

k B T 

hv 


In the Rayleigh-Jeans limit the occupation 


(A5) 


The temperature T oc 7 oc / (not because of a Doppler shift, but from the shock jump conditions), so that 
I v oc / 4 . This result must be divided by f 2 to allow for the fact that the observed radiation was emitted 
at a frequency lower by a factor /, and /„ oc v 2 . Finally, the observed effective radiating area is oc f ~ 2 , 
leading to the result that F v is independent of 7 . 

Each of these derivations uses the shock jump conditions. The final result cannot be obtained from 
Lorentz transformations of the radiation field alone. 
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Fig. 1. - Peak fluxes in four bands for GRB970508. The radio (8.46 GHz) and soft X-ray points are lower 
bounds, indicated by arrows, because the maxima in these bands were not observed. The radio, visible and 
gamma-ray points are not consistent with a single power law. Data are from Kouveliotou et al. (1997) 
(gamma-rays), Costa et al. 1997b (hard X-rays), Piro et al. (1997) (soft X-rays), Galama et al. (1997b) 
and Jaunsen et al. (1997) (V) and Frail et al. (1997b) (radio). The X-ray and gamma-ray points depend 
significantly on assumptions about the spectrum in those bands. 



- 15 - 


gamma-ray 2-30 KeV V Radio 
0.06 1.3 ~ 2 ? 

Table 1: Ratios of peak fluxes, GRB970508 to GRB970228, in four bands. Data from Kouveliotou et al. 
(1997) (extrapolated to match BeppoSAX 40-600 KeV band) and Hurley et al. (1997) (gamma-rays), Costa 
et al. (1997ab) (2-30 KeV X-rays) and from Galama et al. (1997b), Jaunsen et al. (1997) and Groot et al. 
(1997) (V). The entry under V assumes that the characteristic time scale and post-maximum decay rate of 
GRB970228 (whose maximum was not observed) is similar to that of GRB970508; the data directly imply 
a V band ratio of < 3. 


gamma-ray X-ray V 
0.05 ? ? 


1.43 GHz 
> 0.2 


Table 2: Ratios of peak fluxes, GRB970508 to GRB970111, in four bands. Data from Kouveliotou et al. 
(1997) (extrapolated to match BeppoSAX 40-600 KeV band) and Frail et al. (1997a) (gamma-rays) and 
from Frail et al. (1997ab) and Galama et al. (1997a) (radio). A lower bound is given at 1.4 GHz because 
the radio flux of GRB970508 was still rising in the most recent available data. 
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